function [T,n]=bbct(f,~,~,eps)
%f表示被积函数句柄
%a.b表示被积区间[a，b]的输点
% eps 表示精度
%T是用变步长梯形法求得的积分值
%n表示二分区间的次数
h=b-a;
fa=feval(f,a);
fb=feval(f,b);
T1-h*(fa+fb)/2;
T2=T1/2+h*feval(fa-h/2)/2;
n=1;
%按变步长梯形法求积分值；
while abs(T2--T1)>=eps
b=h/2;
Tl=T2;
S=0;
x-alh/2;
while xcb
fx=feval(f,x);
S-S+fx;
x-x-h;
end
T2=T1/2+S*h/2;
n=n+1;
end
T=T2;